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We examine the power spectrum of clusters in the Press-Schechter (PS) 

theory and in N-body simulations to see how the power spectrum of clusters is 






related to the power spectrum of matter density fluctuations in the Universe. 
An analytic model for the power spectrum of clusters for their given number 
density is presented, both for real space and redshift space. We test this model 
■<^|- ' against results from N-body simulations and find that the agreement between 

the analytic theory and the numerical results is good for wavelengths A > QQh~ l 
Mpc. On smaller scales non-linear processes that are not considered in the linear 
PS approximation influence the result. We also use our analytic model to study 
the redshift-space power spectrum of clusters in cold dark matter models with a 
cosmological constant (ACDM) and with a scale-invariant Harrison-Zel'dovich 
initial spectrum of density fluctuations. We find that power spectra of clusters 
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in these models are not consistent with the observed power spectra of the 
APM and Abell-ACO clusters. One possible explanation for the observed 
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power spectra of clusters is an inflationary scenario with a scalar field with the 
potential that has a localized steplike feature. We use the PS theory to examine 
the power spectrum of clusters in this model. 



Subject headings: methods: numerical -galaxies: clustering - galaxies: clusters: 
general - cosmology: theory - large scale structure of universe 



1. INTRODUCTION 

Clusters of galaxies are efficient tracers of the large-scale structure of the universe. A 
strong spatial correlation of clusters of galaxies (Bahcall & Soneira 1983; Klypin & Kopylov 
1983) provided some of the first evidence for the existence of organized structure on large 
scales. To date, much effort has been devoted to determine the correlation function and 
the power spectrum of clusters of galaxies on large spatial scales (e.g. Postman, Huchra 
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& Geller 1992; Peacock & West 1992; Einasto et al. 1993; Dalton et al. 1994; Romer et 
al. 1994; Croft et al. 1997; Einasto et al. 1997a; Retzlaff et al. 1998; Tadros, Efstathiou 
& Dalton 1998). Figure 1 shows the power spectrum of the spatial distribution of the 
Abell-ACO clusters as determined by Einasto et al. (1999), and the power spectrum of the 
APM clusters as found by Tadros, Efstathiou & Dalton (1998). For comparison, we show in 
Figure 1 the power spectrum derived from the distribution of galaxies in the APM survey 
(Baugh & Efstathiou 1993). Observations give the power spectrum of clusters for a given 
number density of clusters, n c \. The number density of the APM clusters and Abell clusters 
is n d ~ 3.4 x 10~ 5 /i 3 Mpc~ 3 and n d ~ 2.5 x 10~ 5 /i 3 Mpc^ 3 , respectively (Dalton et al. 1994; 
Einasto et al. 1997b; Retzlaff et al. 1998). 

As the observational data on the power spectrum of clusters improve, the need for 
precise theoretical predictions becomes increasingly important. The power spectrum of 
clusters for different cosmological models can be calculated using N-body simulations (e.g. 
Retzlaff et al. 1998; Tadros, Efstathiou & Dalton 1998; Suhhonenko & Gramann 1999). 
However, this approach is difficult as it requires simulations with a very large dynamical 
range to identify correctly a sufficient number of clusters. In this situation it would be 
useful to have an analytical approach to describe the origin of the power spectrum of 
clusters of galaxies and to predict it. 

The correlation function of clusters of galaxies has been a subject of many attempts 
at analytical modelling (e.g. Kaiser 1984; Kashlinsky 1987; Cole & Kaiser 1989; Mann, 
Heavens, & Peacock 1993; Mo & White 1996, hereafter MW96; Catelan et al. 1998; 
Kashlinsky 1998). In particular, using the Press-Schechter (PS) formalism to calculate the 
correlation function of clusters in Lagrangian space and mapping from Lagrangian space to 
Eulerian space within the context of the spherical collapse model, MW96 have derived an 
analytical expression for the correlation function of clusters of mass M as 

£ M (r)=& 2 (M)£(r), (1) 

where £(r) is the correlation function of matter density fluctuations. The bias parameter 
b(M) is 

w^+K^fe" 1 )' (2) 

where a 2 (M) is the density dispersion on a mass scale M and 5 t = 1.686. This equation 
has been found to be in good agreement with N-body results by MW96 and by Mo et al. 
(1996). 

The power spectrum and the correlation function of clusters form a Fourier transform 
pair. Therefore, in the linear approximation the power spectrum of clusters of mass M can 
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be expressed as 

P M (k) = b 2 (M)P(k), (3) 

where P(k) is the power spectrum of matter density fluctuations. 

In this paper we develop further the results obtained by MW96 and calculate the 
power spectrum and the correlation function of galaxy clusters for a given number density 
of clusters n c i in Lagrangian space and in Eulerian space (eq. [15,23-24] and [25-27] below). 
Observations provide the distribution of clusters in redshift space, which is distorted due to 
peculiar velocities of clusters. In order to study the power spectrum of clusters in redshift 
space, we use the linear approximation derived by Kaiser (1987). In §3 we examine the 
cluster power spectrum in real space and in redshift space using N-body simulations with 
realistic initial power spectra and show that in both cases the PS theory gives an accurate 
description of the power spectrum of clusters at wavenumbers k < O.lh Mpc^ 1 (or at 
wavelengths A > 60/tT 1 Mpc). In §4 we use our model to study the redshift-space power 
spectrum of clusters in different cold dark matter models. We study inflationary models 
with a scale-invariant initial spectrum of density fluctuations and with a steplike initial 
spectrum derived by Starobinsky (1992). §5 summarizes the main results. 

A Hubble constant of H = lOO/i km s _1 Mpc -1 is used throughout this paper. 



2. THE POWER SPECTRUM OF CLUSTERS IN THE PS 

APPROXIMATION 

An important aspect of the PS model is that, being entirely based on linear theory, 
suitably extrapolated to the collapse time of spherical perturbations, it is, by definition, 
local in Lagrangian space. While this Lagrangian aspect of the theory does not have 
immediate implications for the study of the mass function of clusters of galaxies, it is 
important for the study of their spatial clustering properties. 

Let us assume that the density contrast €m (q) on a mass scale M is a Gaussian random 
field. Here q represents the comoving Lagrangian coordinate. The probability that one 
would find a density contrast between e M and Cm + dtu is 



p(e M ) dt M = ,— — exp 
V27T(T(M 



~M 



2a 2 (M) 



de M . (4) 



The density dispersion <r 2 (M) can be written as 

a 2 (M) =< ^ >= -L | P(k) W 2 (kR) k 2 dk, (5) 
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where W(kR) is the Fourier transform of the window function applied to determine the 
density field. In this paper we will use the top-hat window function. For the top-hat 
window, the mass M is related to the window radius R as M = 47rp&i? 3 /3 (where pi, is the 
mean background density). 



Press & Schechter (1974) suggested that we can determine the probability that a 
cluster of mass M has formed as 



n 



i\i 



p(e M ) de 



M 



- erfc 
2 



St 



y/2a(M)\ ' 



(6) 



where 5 t is a certain threshold value. PS used a value 5 t = 1.686, which is motivated by the 
spherical collapse model. The number density of clusters with the mass between M and 
M + dM is given by (Press & Schechter 1974) 



n{M)dM = 2<^^dM 
y ' M dM 



'2p b S t da{M) 
7i M a 2 (M) dM 



cxp 



s 2 t 



2a 2 (M) 



dM. 



(7) 



Equation (7) involves an additional correction factor 2 to allow all the matter in the 
universe to form bound structures (see Peacock & Heavens 1990; Bond et al. 1991 for a 
more detailed discussion of this correction). Therefore, the number density of clusters of 
mass larger than M can be expressed as 



n d (> M) 



M 



n{M')dM' 



St da(r) 



{2iryi 2 J R a 2 {r) dr 



cxp 



s 2 t 



2a 2 (r) 



dr 



Equation (8) has been frequently used to determine the cluster abundance in different 
cosmological models (e.g. Efstathiou et al. 1988; White, Efstathiou & Frenk 1993; Lacey & 
Cole 1994; Eke, Cole & Frenk 1996; Borgani et al. 1997). 

For a two-dimensional Gaussian distribution the probability density for finding 
simultaneously eu x on scale M\ and em 2 on scale M 2 separated by q is given by 



K e Mi,eM 2 



1 



2na(M 1 )a{M 2 )Jl-p 2 2 



cxp 



2(1 -ph) 



[X-, 



2p 12 x 1 x 2 + x 2 



(9) 



where x± = €mi/o'(Mi), x 2 = €m 2 /°'(M 2 ) and p± 2 = £12(9) / 'cr(Mi)a(M 2 ) is the correlation 
coefficient between €m 1 and €m 2 - The £12(9) is the two-point correlation function of the 
linear density contrast smoothed on the scales Mi and M 2 : 



62(g) =< eM!(qi)eM 2 (q2) >= 



2vr 2 



P(k) W(kR t ) W(kR 2 )j (kq)k 2 dk, 



(10) 
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where q — |qx — q 2 | and Jo(x) is the spherical Bessel function of the order zero. Similarly 
to equation (6), we can write for the probability for two clusters of mass M 1; M 2 to form as 

HmiM 2 = / / p(cm 1 ,CM 2 ) d> e M 1 d€M 2 ■ (H) 

J St J St 

To determine the correlation function of clusters of masses [Mi,Mi + dMi] and [M 2 ;M 2 + dM 2 ] 

we can use an expression 

9 riM 1 M 2 



1 + ^m.mM - du^dn^ - ( 12 ) 

dM\ dM 2 

A similar equation was used by Kashlinky (1987, 1998) to derive the correlation function of 
clusters in Eulerian space and by Catelan et al. (1998) to describe clustering in Lagrangian 
space. 

In order to find the correlation function of clusters of mass larger than M, we start 
from approximations (6-7) and (11-12) and write down the probability to find two clusters 
in small volumes dVi, dV 2 , one cluster being in the mass range [M 1; M\ + dMi], the other 

in [M 2 , M 2 + dM 2 }, as 

d$(Mx, M 2 ) = 4 — ^— % Um ' M2 dM 1 dM 2 dV 1 dV 2 . (13) 

v ' ; M 1 M 2 dM l dM 2 K ' 

We included in equation (13) a correction factor 4 to make it consistent with equation (7) 
at £12(9) = 0. The probability to find two clusters of mass larger than M can be expressed 

as 

fOO fOO 

d$ d (> M)= d$(M u M 2 )dM x dM 2 = n 2 cl [l + £ d (q)] dV x dV 2 , (14) 

Jm Jm 

where n c \ is the number density of clusters of mass larger than M (equation 8). Therefore, 

the correlation function of clusters can be written in the form 

n d Jm Jm M 1 M 2 oMiOM 2 

Equations (5,8-11) and (15) determine the correlation function of clusters for a given n c i in 
the PS theory of gravitational clustering. 

Let us consider the correlation function of clusters in the linear approximation, i.e. 
neglecting the terms involving ^\ 2 {q) in equation (9). In the linear approximation, the 
probability density 



p(cmi , £m 2 



27ra(M 1 )a(M 2 ) 



, £l2 (^M^Ma 



a 2 (M!)a 2 (M 2 ) 



exp 
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2a 2 (Mi 
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€ M 2 



2a 2 (M 2 



(16) 



-6- 



Substituting approximation (16) into equation (11) we find that 

1 62(g) 



n 



MiM 2 



^■mJ^m 2 + 



and 



d 2 U 



M\M 2 



dU Ml dU 



2tt a(M 1 )a(M 2) 
62(g) 



exp 



s? 



2a 2 (Mi 



exp 



^ 



2a 2 (M 2 



(17) 



iii 2 



dM l dM 2 dM x dM 2 



1 + 



^ 



5} 



o 2 (M 2 



5 2 \a 2 (M 1 ) 

For separations q much larger that the smoothing lengths, q » Ri and q » R 2 , the 
influence of the window functions on the correlations is negligible, and the correlation 
function £,\ 2 (q) — C,(q) (where £(g) is the linear mass autocorrelation function). Therefore, 
in the linear approximation the correlation function of clusters of masses [Mi\Mi + dMi] 
and [M 2 ;M 2 + dM 2 ] can be written as 



ZM 1 MM = b L (M 1 )b\M 2 )Z(q), 
where the Lagrangian bias factor b L (M) is 

1 I s 2 



b L (M) 



S t 



1 



(19) 



(20) 



a 2 (M) 

MW96 derived this equation to describe clustering in Lagrangian space starting from the 
behaviour of the conditional Lagrangian mass function derived by Bond et al. (1991). By 
using the spherical collapse model, MW96 found that in Eulerian space the correlation 
function of clusters of masses [MijMi + dM\\ and [M 2 ;M 2 + dM 2 ] can be written as 



tM 1 M 2 (r) = b(M 1 )b(M 2 )£(r), 



(21) 



were £(r) is the correlation function of matter density fluctuations in Eulerian space and 
the Eulerian bias factor b(M) is 



b(M) = l + b L (M) = 1 + 



tf 



St \a 2 (M) 



1 



(22) 



Catelan et al. (1998) used a more general method for evolving the spatial distribution of 
clusters and found a similar result for the linear approximation. The shift by 1 of the linear 
bias factor, caused here by the transformation from the Lagrangian to the Eulerian world, 
comes from mass conservation in Eulerian space. 

Let us consider the correlation function of clusters for a given number density in the 
linear PS approximation. Substituting approximation (18) into equation (15) we find that 
in the linear approximation the Lagrangian correlation function of clusters for a given n c ; 
can be expressed as 

UQ) = (&5) 2 *(?) , (23) 
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where the bias parameter b d can be written in the form 



n cl JM 



b L {M')n{M')dM' 



1 da{r) 
(2n) 3 / 2 n c iJR cr 2 (r) dr 



% 



a 2 (r) 



cxp 



% 



2a 2 (r) 



dr 
■(24) 



Similarly to equations (21-22) the Eulerian correlation function of clusters for a given 
number density n c i can be written as 



Ur) = b 2 d d{r) 



(25) 



where the bias parameter b c i is 

3 



1 + 6, 



1- 



da{r) 



(2n) 3 / 2 n c iJR <j 2 (r) dr 



a 2 (r) 



-1 



cxp 



% 



2a 2 (r) 



dr 



(26) 



The power spectrum of clusters for a given number density in the linear approximation can 
be expressed as 

P d (k) = b d P(k). (27) 

The cluster bias parameter b c i depends on the minimal mass M (or the window radius R) 
of clusters and on the power spectrum of density fluctuations, P(k), which determines the 
function <r(r). For a fixed P(k) and n d , the minimal mass M (or scale R) can be determined 
by inverting equation (8). In this approach the power spectrum (or the correlation function) 
of clusters does not depend on the mean background density pb (or the density parameter 
fi ). 

Observations provide the distribution of clusters in redshift space, which is distorted 
due to peculiar velocities of clusters. On large scales, where linear theory applies, the power 
spectrum of matter density fluctuations in redshift space is given by (Kaiser 1987): 



P s (k) 



1 



2/(fi ) , PWj 



5 



P(k), 



(28) 



where /(^o) ~ ^o' 6 is the linear velocity growth factor. In the linear approximation (27), 
relation (28) takes the form 



P d (k) 



x + 2/(^0) + p(n ) 



3b r 



5^ 



b d P(k). 



(29) 



Equation (29) determines the power spectrum of clusters for a given n c i in redshift space. 



3. NUMERICAL RESULTS 

For testing the solutions (26-27) and (29), we ran N-body simulations, using a 
particle-mesh code described by Gramann (1988). We investigated the evolution of 256 3 
particles on a 256 3 grid, with $1 = 1. The comoving box size was L = 384/i" 1 Mpc. The 
initial density field was taken to be Gaussian. 

We examined the distribution of clusters in two cosmological models which start from 
the observed power spectra of the distribution of galaxies and clusters of galaxies. In the 
model (1), the initial linear power spectrum of density fluctuations was chosen in the form 
P(k) oc k~ 2 at wavelengths A < 120/i _1 Mpc. In the model (2), we assumed that the initial 
power spectrum contains a primordial feature at the wavelengths A ~ 30 — 60/i -1 Mpc that 
correspond to the scale of superclusters of galaxies. Suhhonenko & Gramann (1998) studied 
the mass function, peculiar velocities, the power spectrum and the correlation function of 
clusters in both models for different values of the density parameter $lo and as (the rms 
fluctuation on the 8/i _1 Mpc scale). The results were compared with observations. They 
found that in many aspects the initial power spectrum of density fluctuations in the model 
(2) fits the observed data better than the simple power law model (1). 

Clusters were selected in the simulations as maxima of the density field that was 
determined on a 256 3 grid using the CIC-scheme. To determine peculiar velocities of 
clusters, we determined the peculiar velocity field on a 256 3 grid using the CIC-scheme 
and found the peculiar velocties at the grid points were clusters had been identified. 
The clusters were then ranked according to their density and we selected N c i = (L/d c i) 3 
highest ranked clusters to produce cluster catalogs with a mean intercluster separation 
d c i = 30 — 40/i _1 Mpc. For comparison, the mean separation of the observed APM and Abell 
clusters is d c i ~ 31/i -1 and d c i ~ 34/i _1 Mpc, respectively (Dalton et al. 1994; Retzlaff et al. 
1998). 

It is difficult to follow the evolution of rich clusters by using N-body simulations, as 
it requires simulations with a very large dynamical range to identify correctly a sufficient 
number of clusters. We determined clusters as maxima of the density field smoothed on 
the scale R ~ 1.5/i _1 Mpc. This method of identifying clusters is not really identical to 
that what the PS theory is predicting - the location of collapsed, virialized objects - and 
this could affect the degree of agreement between the analytic theory and the numerical 
results. However, in order to increase the resolution of simulations we must increase the 
number of test particles and grid points, or we have to follow the evolution of clusters in a 
smaller box. While the first possibility is technically difficult, in the latter case the number 
of rich clusters becomes too small to get statistically reliable results. Taking into account 
the requirements on the number of clusters and on the resolution together with the fact 
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of fixed computer recources we decided to use a box size L = 384/i 1 Mpc and a grid size 
R g = 1.5/i- 1 Mpc. 

Figure 2 shows the power spectrum of clusters with a mean separation d c \ = 30/i _1 
Mpc in our models. We studied also the power spectrum of clusters with a mean separation 
dd = 35/i~ 1 Mpc and d c \ = 40/i" 1 Mpc, and found similar results. To calculate the power 
spectrum of clusters in the simulations we determined the density field of clusters on a 128 3 
grid using the CIC-scheme and calculated its Fourier components, subtracting the shot 
noise term. We determined also the Poisson error bars for the power spectrum. In a shell of 
/c-space containing m modes, the Poisson error can be estimated as AP c i(k) = 2 1 / 2 vrr 1 l 2 d'^. l 
(see e.g. Peacock & Nicholson (1991) for a more detailed analysis of the Poisson errors in 
the power spectrum). 

Figure 2a shows the results for the model (1) and Figure 2b for the model (2). For the 
model (1), the clusters were determined in the simulation for a$ = 0.5 and a§ = 0.8. For 
the model (2), these parameters were a§ = 0.5 and <r 8 = 0.84. For comparison, we examined 
the power spectrum of clusters in the linear PS approximation (27). We inverted equation 
(8) to determine the minimal radius R for a given n c i (or d c i) and after that we calculated 
the bias factor b c i using equation (26). For a given og and n c i, the bias parameter b c \ in 
the model (1) is similar to that in the model (2). For <jg = 0.5, we find that b c \ = 4.3 and 
4.9 for separations d c i = 30 and 40/i~ 1 Mpc, respectively. For erg = 0.8, b c i = 2.9 and 3.3, 
respectively. 

Now we can compare the numerical results and the PS theory predictions. First, 
consider the power spectrum of clusters at wavenumbers k = 0.04 — O.lh Mpc -1 (or 
at wavelengths A = 60 — 160/i^ 1 Mpc). On larger scales, the number of modes in the 
simulations is too small to get statistically reliable results. Figure 2 demonstrates that at 
wavenumbers k = 0.04 — O.lh Mpc -1 , the power spectrum of clusters in the simulations is 
linearly enhanced with respect to the power spectrum of the matter distribution in both 
models studied. For <jg = 0.5, the agreement between the results of the N-body simulations 
and the PS theory predictions is very good. The mean deviation between the numerical 
results and the theoretical predictions is about 4% and 7% in the model (1) and in the 
model (2), respectively. For as ~ 0.8, we find that the power spectrum of clusters in N-body 
simulations is somewhat lower than predicted by approximation (27). The power spectrum 
of clusters is about 80% and 70% of the linear theory predictions in the model (1) and 
model (2), respectively. The PS approximation predicts that during the evolution between 
the as = 0.5 and a 8 = 0.8, P c i(k) increases slightly. In simulations we find that during 
the evolution the power spectrum of clusters decreases. This effect is probably caused by 
merging of very rich clusters. Further study (e.g. numerical simulations with a higher 
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dynamical range) is needed to determine whether this is a real effect in our models, or a 
numerical effect due to the limited dynamical range of the N-body simulations. 

Let us now consider the power spectrum of clusters at wavenumbers k = 0.1 — 0.2ft, 
Mpc -1 (A = 30 — 60/i -1 Mpc). In the model (1), the numerical results on these scales are 
in good agreement with the linear approximation (27). But in the model (2), the power 
spectrum of clusters in the simulation is significantly smaller than that predicted by the 
linear theory. For a§ = 0.84, we find that the power spectrum of clusters is only 30% 
of the linear theory predictions. Therefore, in the model (2), we cannot use the linear 
approximation to study the power spectrum of clusters at wavenumbers k > O.lh Mpc -1 . 
The power spectrum of clusters on these scales is probably determined by non-linear 
processes in superclusters of galaxies. 

Figure 3 shows the redshift-space power spectrum of clusters in the simulations and 
in approximation (29). For clarity, we did not plot error bars in this Figure. In order 
to study peculiar velocities of galaxy clusters and their distribution in redshift space in 
models with different Qq, we determined peculiar velocities of clusters in the simulations 
with Q = 1 and assumed that peculiar velocities of clusters are proportional to the linear 
growth factor f(£lo). Let us consider the redshift-space power spectrum of clusters at 
wavenumbers k = 0.04 — O.lh Mpc -1 , where we expect that the linear approximation (29) 
applies. Numerical results show that on these scales the power spectrum of clusters in 
redshift space is linearly enhanced with respect to the power spectrum of clusters in real 
space. In Figure 3 we show the power spectrum of clusters in redshift space for o" 8 = 0.5 and 
Q = 1. In this case, the agreement between the numerical results and approximation (29) is 
very good. The mean deviation between the numerical results and theoretical predictions is 
about 5% and 7%, in the model (1) and model (2), respectively. For a 8 rs 0.8 we find that, 
similarly to real space, the redshift-space power spectrum of clusters in the simulations is 
somewhat lower than that predicted by the linear theory. For Q = 1, the power spectrum of 
clusters is about 80% and 75% of the linear theory predictions in the model (1) and model 
(2), respectively. For Q Q = 0.2, we found similar results. Within the uncertainties due to 
various numerical effects in the simulations, the agreement between the analytic theory and 
the numerical results is good. The degree of agreement is similar in real and in redshift 
space. 



4. MODELS WITH COLD DARK MATTER 

Now we use approximation (29) to analyze the redshift-space power spectrum of 
clusters in different cold dark matter models. We examine flat cosmological models with 
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the density parameter Q = 0.3 — 0.4, the baryonic density Qb — 0.015/i -2 and the 
normalized Hubble constant h = 0.5 — 0.7. These parameters are in agreement with recent 
nucleosynthesis results (Tytler at al. 1996), with measurements of the density parameter 
(e.g. Dekel, Burstein & White 1996; Bahcall & Fan 1998) and with measurements of the 
Hubble constant using various distance indicators (e.g. Tammann 1998). To restore the 
spatial flatness in the low-density models, we assume a contribution from a cosmological 
constant: Q\ = 1 — Sl - 

Figure 4 shows the redshift-space power spectrum of clusters in the ACDM models with 
a scale-invariant Harrison-Zel'dovich initial spectrum of density fluctuations (Pi n (k) oc k). 
We have used the transfer function derived by Bardeen et al. (1986) and Sugiyama (1995), 
and the COBE normalization derived by Bunn & White (1997). Figure 4a demonstrates the 
power spectrum of the model clusters and the Abell-ACO clusters with a mean separation 
d c i = 34/i -1 Mpc. We show the power spectrum of the Abell-ACO clusters as found by 
Einasto et al. (1999). This spectrum represents the weighted mean of the power spectra 
determined by Einasto et al. (1997a) and Retzlaff et al. (1998). Einasto et al. (1997a) 
determined the power spectrum of the Abell-ACO clusters from the correlation function of 
clusters, while Retzlaff et al. (1998) estimated the power spectrum directly (see Einasto et 
al. 1999 for details). The power spectrum of the distribution of the Abell-ACO clusters 
peaks at the wavenumber k = 0.052/i Mpc -1 (or at the wavelength A = 120/i -1 Mpc). 
For k > 0.052/i Mpc -1 , the cluster power spectrum is well approximated by a power law, 
P{k) oc k n , with n rs —1.9. A similar peak in the one-dimensional power spectrum of a deep 
pencil-beam survey was detected by Broadhurst et al. (1990) and in the two-dimensional 
power spectrum of the Las Campanas redshift survey by Landy et al. (1996). Figure 4b 
shows the power spectrum of the model clusters with a mean separation d c i = 31/i -1 Mpc. 
For comparison, we show the power spectrum of the observed APM clusters determined 
by Tadros, Efstathiou & Dalton (1998). They analyzed the redshift survey of 364 clusters 
described by Dalton et al. (1994). The mean intercluster separation of APM clusters is 
d d ~ 31k- 1 Mpc (Dalton et al. 1994). 

Figure 4 shows that the power spectrum of clusters predicted in the ACDM models 
is not consistent with the observed spectra of clusters. We studied a x 2 probability at 
wavenumbers k = 0.03 — O.lh Mpc -1 , where we expect that the linear approximation (29) 
applies. For the models presented in Figure 4, the probability to fit the observed power 
spectra of the Abell and APM clusters is less than 5 x 10 -2 and 10 -6 , respectively. We 
also studied a \ 2 probability for the models, where the amplitude of the power spectrum of 
clusters is about 70% and 80% of the linear theory predictions and found similar results. 
The power spectrum of clusters in the ACDM models with a scale- invariant initial spectrum 
is not consistent with the observed spectra of the APM and Abell-ACO clusters (see also 
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e.g. Tadros, Efstathiou & Dalton 1998; Einasto et al. 1997a). 

The power spectrum of density fluctuations in the universe depends on the physical 
processes in the early universe. The peak in the power spectrum of clusters at wavelength 
A ~ 120/i" 1 Mpc may be generated during the era of radiation domination or earlier. 
Baryonic acoustic oscillations in adiabatic models may explain the observed power spectrum 
only if currently favored determinations of cosmological parameters are in substantial error 
(e.g. the density parameter Q < 0.2h) (Eisenstein et al. 1998). One possible explanation 
for the observed power spectra of clusters is an inflationary model with a scalar field whose 
potential V(tp) has a local steplike feature in the first derivative. This feature can be 
produced by fast phase transition in physical field other than an inflaton field. An exact 
analytical expression for the scalar (density) perturbations generated in this inflationary 
model was found by Starobinsky (1992) (see also Lesgourgues, Polarski & Starobinsky 
1998). The initial power spectrum of density fluctuations in this model can be expressed as 

Ptn(k) K »WM ( 3o) 

where function S(k/ko,p) can be written in the form 



S(y,p) = 1 

y 



2 
/i(y)sin2y + -cos2y 

y 



9p?/ 2 (y) 



2 
h(y) + fi (y) cos 2y sin 2y 

y 



2y 2 

"(31) 

Here, function fi(y) = 1 — y~ 2 , /2(y) = 1 + y~ 2 and p\ = p — 1. The initial power spectrum 
in this model depends on two parameters ko and p. The parameter ko determines the 
location of the step and the parameter p - the shape of the initial spectrum. For p = 1, 
we recover the scale-invariant Harrison-Zel'dovich spectrum (S(y, 1) = 1). At present, the 
initial spectrum (30-31) is probably the only example of a initial power spectrum with the 
desired properties, for which a closed analytical form exists. 

Figure 5 shows the power spectrum of clusters in the ACDM models with a steplike 
initial power spectrum. As in Figure 4, we have used the transfer function derived by 
Bardeen et al. (1986) and Sugiyama (1995), and the COBE normalization derived by Bunn 
& White (1997). For the models with p < 1 and p > 1, the step parameter was chosen 
to be ko = O.OI6/1 Mpc -1 and ko = 0.03/i Mpc" 1 , respectively. In this case, in the models 
with p < 1, the power spectrum has a well-defined maximum at the wavenumber k ~ 0.05/i 
Mpc" 1 and a second maximum at k ~ O.l/i Mpc -1 . In the models with p > 1, the picture is 
inverted. The power spectrum has a flat upper plateau at wavenumbers k < 0.05/i Mpc" 1 , 
a sharp decrease on smaller scales (k = 0.05 — O.lh Mpc" 1 ) and a secondary maximum at 
k ~ 0.15/i Mpc" 1 . At wavenumbers k > 0.05h Mpc -1 , the power spectrum in the ACDM 
models with p ~ 1.3 — 1.4 is similar to the power spectrum in the numerical model (2), that 
we studied in §3 (see Figure 2b and Figure 3b). 
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Figure 5a shows the power spectrum of the model clusters and the Abell clusters with 
a mean separation d c i = 34/i _1 Mpc. We have studied inflationary models with parameter 
p = 0.6 — 0.8. Figure 5a shows that the shape of the cluster power spectrum in these 
models is in good agreement with the observed power spectrum of the Abell- ACO clusters. 
However, the amplitude of the observed spectrum of clusters is about 70% lower than 
predicted in these models by using approximation (29). We found a similar effect in §3 by 
comparing the analytic and numerical results. It is possible that the linear approximation 
(29) somewhat overestimates the power spectrum of clusters due to dynamical effects that 
are not taken into account in this approximation. Therefore, the power spectrum of clusters 
can be expressed as 



pzm 



i | 2/(n ) | p(n ) 



3b d 5bl 



KiP(k), (32) 



where the factor F = 0.7 — 1.0. For F = 0.7, the inflationary models studied in Figure 5a, 
are consistent with the power spectrum of the Abell- ACO clusters at a confidence level of 
> 90% (based on a x 2 test at k — 0.03 — O.lh Mpc -1 ). Figure 5b demonstrates the power 
spectrum of the model clusters and the APM clusters with a mean separation da = 31/i 1 
Mpc. We have studied the ACDM model with p = 1.25. For F = 0.7, the cluster power 
spectrum in this model is consistent with the observed power spectrum of the APM clusters 
at a confidence level of > 50%. (For F = 0.65, at a confidence level of > 90%). 

Available data are insufficient to rule out any of the models studied in Figure 5. Note 
that the power spectra of clusters predicted in the models with p < 1 and p > 1 are rather 
different on large scales, where k < 0.05/i Mpc -1 . Therefore, accurate measurements of the 
power spectrum of clusters and galaxies on these scales can serve as a discriminating test 
for these interesting models. 



5. SUMMARY AND CONCLUSIONS 

Because of the low density of rich clusters of galaxies and their enhanced clustering 
strength compared to galaxies, rich clusters provide a powerful probe of large-scale structure 
in the Universe. In this paper, we have examined the power spectrum of clusters in the 
Press-Schechter (PS) theory and in N-body simulations to see how the power spectrum of 
clusters is related to the power spectrum of matter density fluctuations in the Universe. We 
have derived an analytical expression to determine the correlation function of clusters for 
a given number density of clusters in Lagrangian space and have examined this expression 
in the linear approximation. In order to study the power spectrum of clusters in redshift 
space, we used the approximation derived by Kaiser (1987). 
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For testing our analytic results, we used N-body simulations with realistic initial power 
spectra of density fluctuations. We investigated the power spectrum of clusters with a mean 
separation d c i = 30 — 40/i -1 Mpc. The numerical results showed that we can use the linear 
PS approximation to predict the power spectrum of clusters in real space and in redshift 
space at wavenumbers k < O.lh^ 1 Mpc. On smaller scales non-linear processes that are not 
taken into account in the linear PS approximation influence the results. 

We also used the PS approximation to analyze the redshift-space power spectrum of 
clusters in ACDM models. We investigated inflationary models with a scale-invariant initial 
spectrum of density fluctuations and with a steplike initial spectrum derived by Starobinsky 
(1992). The results were compared with observations. We found that the power spectrum 
of clusters in the ACDM models with a scale-invariant initial spectrum is not consistent 
with the observed spectra of the APM and Abell clusters. We investigated also inflationary 
models with a steplike initial power spectrum (30-31) that depends on two parameters ko 
and p. The parameter ko determines the location of the step and the parameter p - the 
shape of the initial spectrum. For the models with p < 1 and p > 1, the step parameter 
was chosen to be k = 0.016/i Mpc" 1 and k = 0.03h Mpc -1 , respectively. We found that 
the power spectrum of clusters in these ACDM models is in good agreement with the 
observed power spectrum of the Abell-ACO clusters, if the initial parameter p is in the 
range p = 0.6 — 0.8. To describe the power spectrum of the APM clusters, we can use the 
ACDM model with parameter p = 1.25. 

We thank J. Einasto, A. Kashlinsky, S. Matarrese, H. Mo and E. Saar for useful 
discussions. This work has been supported by the ESF grant 97-2645. 



-15- 



REFERENCES 

Bahcall, N. A., & Fan, X. 1998, Proc. Nat. Acad. Sci. USA, 95, 5956 

Bahcall, N. A., & Soneira, R. M. 1983, ApJ, 270, 20 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 

Baugh, C, & Efstathiou, G. 1993, MNRAS, 265, 145 

Bond, J. R., Cole, S., Efstathiou, G., & Kaiser, N. 1991, ApJ, 379, 440 

Borgani, S. et al. 1997, NewA, 1, 321 

Broadhurst, T.J., Ellis, R.S., Koo, D.S., & Szalay, A.S. 1990, Nature, 343, 726 

Bunn, E.F., & White, M. 1997, ApJ, 480, 6 

Dalton, G. B., Croft, R. A. C, Efstathiou, G., Sutherland, W. J., Maddox, S. J., & Davis, 
M. 1994, MNRAS, 271, L47 

Dekel, A., Burstein, D., & White, S. D. M. 1996, in Critical Dialogues in Cosmology, ed. 
Neil Turok, (Singapore: World Scientific) 

Catelan, P., Lucchin, F., Matarrese, S., & Porciani, C. 1998, MNRAS, 297, 692 

Cole, S., & Kaiser, N., 1989, MNRAS, 237, 1127 

Croft, R. A. C, Dalton, G. B., Efstathiou, G., Sutherland, W. J., & Maddox, S. J. 1997, 
MNRAS, 291, 305 

Efstathiou, G., Frenk, C.S., White, S.D.M., & Davis, M. 1988, MNRAS, 235, 715 

Einasto, J., Gramann, M., Saar E., & Tago, E. 1993, MNRAS, 260, 705 

Einasto, J. et al. 1997a, Nature, 385, 139 

Einasto, M., Tago, E., Jaaniste, J., Einasto, J., & Andernach, H. 1997b, A&AS, 123, 119 

Einasto, J. et al. 1999, ApJ, in press (astro-ph: 9812247) 

Eisenstein, D.J., Hu, W., Silk, J., Szalay, A.S. 1998, ApJ, 494, LI 

Eke, V. R., Cole, S., & Frenk, C. S. 1996, MNRAS, 282, 263 

Gramann, M. 1988, MNRAS, 234, 569 

Kaiser, N. 1984, ApJ, 284, L9 

Kaiser, N. 1987, MNRAS, 227, 1 

Kashlinsky, A. 1987, ApJ, 317, 19 

Kashlinsky, A. 1998, ApJ, 492, 1 

Klypin, A. A., & Kopylov, A. I. 1983, Soviet. Astron. Lett., 9, 41 



-16- 



Lacey, C, & Cole, S. 1994, MNRAS, 271, 676 

Landy, S.D., Shectman, S.A., Lin, H., Kirshner, R.P., Oemler, A., Tucker, D., & Schechter, 
P.L. 1996, ApJ, 456, LI 

Lesgourgues, J., Polarski, D., Starobinsky, A.A., 1998, MNRAS, 297, 769 

Mann, R.G., Heavens, A. F., & Peacock, J. A., 1993, MNRAS, 263, 798 

Mo, H. J., Jing, Y. P., & White, S. D. M., 1996, MNRAS, 282, 1096 

Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 347 (MW96) 

Peacock, J. A., & Heavens, A. F. 1990, MNRAS, 243, 133 

Peacock, J. A., & Nicholson, D., 1991, MNRAS, 253, 307 

Peacock, J. A., & West, M. 1992, MNRAS, 259, 494 

Postman, M., Huchra, J., & Geller, M. 1992, ApJ, 384, 404 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Retzlaff, J., Borgani, S., Gottlober, S., Klypin, A., & Miiller, V. 1998, NewA, 3, 631 

Romer, A. K., Collins, C. A., Cruddace R. C, MacGillivray, H., Ebeling, H., Bohringer, H. 
1994, Nature, 372, 75 

Starobinsky, A. A., 1992, JETP Lett., 55, 477 

Sugiyama N. 1995, ApJS, 100, 281 

Suhhonenko, I. & Gramann, M. 1999, MNRAS, 303, 77 

Tadros, H., Efstathiou, G., & Dalton, G. B. 1998, MNRAS, 296, 995 

Tammann, C, 1998, in General Relativity, 8th Marcel Crossmann Symposium, ed. T. 
Piran, (Singapore World Scientific) 

Tytler, D., Fan, X. M., & Buries, S. 1996, Nature, 381, 207 

White, S. D. M., Efstathiou, G., & Frenk, C. S. 1993, MNRAS, 262, 1023 



This preprint was prepared with the AAS IATj^X macros v4.0. 



-17- 

FIGURE CAPTIONS 

Fig. 1. — The power spectrum of the distribution of clusters of galaxies. Open circles 
show the power spectrum of the distribution of the Abell-ACO clusters and open triangles 
represent the power spectrum of the APM clusters. Filled circles show the power spectrum 
of the galaxy distribution in the APM survey. 

Fig. 2. — The power spectrum of the distribution of clusters in real space. The heavy lines 
show the power spectra of clusters determined by approximation (27) and the light lines 
the corresponding linear power spectra of matter density fluctuations. Panel (a) shows the 
power spectrum in the model (1) for <t 8 = 0.5 (dot-dashed lines) and for a 8 = 0.8 (solid 
lines). Open circles and triangles show the power spectra of clusters in N-body simulations 
for ag = 0.5 and a 8 = 0.8, respectively. Panel (b) shows the power spectrum in the model 
(2) for 0- 8 = 0.5 (dot-dashed lines) and for a 8 = 0.84 (solid lines). Open circles and triangles 
show the numerical results for a 8 = 0.5 and erg = 0.84, respectively. Error bars in (a) and (b) 
denote Poisson errors. The power spectrum is shown for the clusters with a mean separation 
d cl = 30/i- 1 Mpc. 

Fig. 3. — The power spectrum of the distribution of clusters in redshift space. Panel (a) shows 
the results in the model (1) and panel (b) in the model (2). The heavy lines show the power 
spectra of clusters determined by approximation (29) and the light lines the corresponding 
linear power spectra of matter fluctuations. Symbols describe the power spectra of clusters 
in the simulations. In the model (1), we studied clusters for a 8 = 0.5, Q = 1 (dot-dashed 
lines, open circles), for er 8 = 0.8, ffo = 0.2 (dotted lines, crosses) and for a 8 = 0.8, Q = 1 
(solid lines, open triangles). In the model (2), clusters were studied for a 8 = 0.5, Q = 1 (dot- 
dashed lines, open circles), for a 8 = 0.84, Q = 0.2 (dotted lines, crosses) and for a 8 = 0.84, 
Q = 1 (solid lines, open triangles). The power spectrum is shown for the clusters with a 
mean separation d c i = 30/i -1 Mpc. 

Fig. 4. — The redshift-space power spectrum of clusters in the ACDM models with Qq = 0.4, 
h = 0.7 (solid lines), Qq = 0.3, h = 0.7 (dot-dashed lines), Qq = 0.4, h = 0.5 (dashed 
lines) and Qq = 0.3, h = 0.5 (dotted lines). The initial spectrum is assumed to be scale- 
invariant (Pi n (k) oc k). The heavy lines show the power spectra of clusters determined by 
approximation (29) and the light lines the corresponding linear power spectra of matter 
fluctuations, (a) The power spectrum of the model clusters and the Abell-ACO clusters 
(open circles) with a mean separation d c i = 34/i _1 Mpc. (b) The power spectrum of the 
model clusters and the APM clusters (open triangles) with a separation d c i = 31/i _1 Mpc. 

Fig. 5. — The redshift-space power spectrum of clusters in the ACDM models with a steplike 
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initial power spectrum. The heavy lines show the power spectra of clusters determined by 
approximation (32) and the light lines the corresponding linear power spectra of matter 
fluctuations, (a) The power spectrum of the model clusters and the Abell-ACO clusters 
(open circles) with a separation d c \ = 34/i _1 Mpc. We have studied the models with Qq = 0.3, 
h = 0.6, p = 0.8 (solid lines) and Q = 0.3, h = 0.5, p = 0.6 (dashed lines). The cluster 
power spectrum is determined for F = 0.7. The dot-dashed line shows the cluster power 
spectrum in the first model for F = 1.0. (b) The power spectrum of the model clusters 
and the APM clusters (open triangles) with a mean separation d c i = 31/i -1 Mpc. We have 
studied the model with Q = 0.4, h = 0.6 and p = 1.25. The power spectum of clusters is 
determined for F = 0.7 (solid line) and F — 1.0 (dot-dashed line). 
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